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Abstract 

A novel crystalline structure of hybrid monolayer hexagonal boron nitride (BN) and graphene is predicted by means of 
the first-principles calculations. This material can be derived via boron or nitrogen atoms which are substituted by 
carbon atoms evenly in the graphitic BN with vacancies. The corresponding structure is constructed from a BN 
hexagonal ring linking an additional carbon atom. The unit cell is composed of seven atoms, three of which are boron 
atoms, three are nitrogen atoms, and one is a carbon atom. It shows a similar space structure as graphene, which is 
thus coined as g-BsNsC. Two stable topological types associated with the carbon bond formation, i.e., C-N or C-B 
bonds, are identified. Interestingly, distinct ground states of each type, depending on C-N or C-B bonds, and 
electronic bandgap as well as magnetic properties within this material have been studied systematically. Our work 
demonstrates a practical and efficient access to electronic properties of two-dimensional nanostructures, providing 
an approach to tackling open fundamental questions in bandgap-engineered devices and spintronics. 
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Background 

Two-dimensional (2D) nanomaterials, such as graphene 
and monolayer hexagonal boron nitride (h-BN), are 
expected to play a key role in future nanotechnology 
as well as to provide potential applications in next- 
generation electronics. Recently, novel hybrid structures 
consisting of a patchwork of BN and C nanodomains 
(BNC) were synthesized through the use of a thermal cat- 
alytic chemical vapor deposition method [1]. This finding 
immediately has attracted a great deal of research interest 
[2-4], given that it demonstrates a hitherto efficient route 
to tune the bandgaps of these 2D materials. 

It is well known that the perfect hexagonal and planar 
structure of BNC largely depends on the good match- 
ing between BN and C domains. However, it is indeed an 
outstanding challenge as BN and C phases are naturally 
immiscible in 2D [1]. This explains why Ci et al. [1] could 
observe some wrinkles in the atomic force microscopy 
image. Such mutual contradiction mainly originates from 
the domain boundary effect and the staggered potentials 
of B and N atoms in BNC, which doubtlessly affects their 
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continuous tunable electronic energy gaps. It has been 
confirmed in the related theoretical calculations [5-9] 
where the bandgaps show a strong oscillation feature. 

The X-ray photoelectron spectroscopy of BNC mea- 
sured in the work of Ci et al. [1] shows two additional types 
of C-bonding configurations, which correspond to the C- 
B and C-N bonds with the bonding energies of around 
188.4 and 398.1 eV, respectively. This feature means that 
two inequivalent C-bonding types, i.e., C-B or C-N bonds, 
must be present in the boundaries of the hybridized BN 
and C domains, which have a significant effect on deter- 
mining the local structures and subsequently vary the 
electronic properties in this system. For example, the 
transport channels show a robust characteristic gap when 
the topological index changes the sign of the valley Hall 
effect [3]. In addition, both Raman D band at 1,360 cm -1 
and D' band at 1,620 cm -1 are also observed in BNC [1], 
which were attributed to the lattice disorder or the finite 
crystal size. This lattice disorder effect might directly 
introduce vacancies to this 2D hexagonal system [10-12], 
which is also true in BNC as shown in Figure two (a) of [1]. 

For the planar BNC structure, although the larger 
domains would be preferred to decrease the total 
domain interfacial energy, the randomly distributed 
hybrid domains and the immiscible phases, as well as 
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the induced vacancies, must result in various complex 
structures of BNC [13-16]. An immediate consequence 
is the largely inaccessible synthesis of expected BNC in 
experiments [17], which next hinders realization of 
bandgap-engineered applications in actual devices. In this 
endeavour, exploring the structures and electronic prop- 
erties associated with C bond formation in BNC which 
contains C-N or C-B bonds is an interesting topic that 
must be addressed before widespread synthetic applica- 
tions. Thus, a simple model of BNC where the C bonds 
play a crucial role must be considered again. More impor- 
tantly, a deep theoretical understanding, which originally 
was concealed behind the complex hybridized structures, 
is imperative. 

Here, we report that such a simple model of BNC may be 
just the graphitic B3N3C (g-B3N3C), which is a perfect 2D 
monolayer graphite-like structure as shown in Figure 1. 
As mentioned above, the C atom can bond to B and N 
atoms. Therefore, two topological types of g-BsNsC are 
easily deduced. One is the a-g-B^NsC which is related to 
the the higher bonding energy of the C-N bond, while the 
other one is the fi-g-B^NsC which is constructed based on 
the lower bonding energy of the C-B bond. It can be seen 
that such a material is essentially a C-doped graphitic BN 
(g-BN) with vacancies [11]. The substitution of the N/B 
atom with a C atom in g-BN with the B/N vacancy will 
yield the a//?-g-E>3N3C structure. The interactions among 
the C atoms and/or vacancies as well as the C-bonding 
types (C-B and C-N bonds) ing-BsNsC significantly alter 
its electronic properties. To explore this effect, standard 
density functional theory with different functional (see 




4.55 4.60 4.65 4.70 4.75 5.20 5.25 5.30 5.35 

Lattice constant a (A) 



Figure 1 The total energy per atom as a function of lattice 
constant a for g-B 3 N 3 C. Yellow, green, and gray balls represent B, N, 
and C atoms, respectively. Their respective 2x2 supercells are also 
given nearby. 



the following discussion) calculations has been carried 
out for this predicted material. Remarkably, such material 
displays two distinct electronic structure properties: a-g- 
B3N3C is a semiconductor, while fi-g-B^NsC behaves like 
a metal and leads to a magnetic ground state. 

This paper is arranged as follows: In the second section, 
we present the computational method used in this work, 
followed by the electronic bandgap in a-g-BsNsC and 
magnetism in p-g-B^NsC in the third section. We then 
conclude this paper in the fourth section. 

Methods 

For structural optimization, we employed density func- 
tional theory with the generalized gradient approxima- 
tion (GGA) of Perdew-Burke-Ernzerhof (PBE) [18] for 
the exchange-correlation (XC) potential within the pro- 
jector augmented wave method as implemented in VASP 
[19,20]. An all-electron description, the projector aug- 
mented wave method, is used to describe the electron-ion 
interaction. The cutoff energy for plane waves is set to 
be 500 eV, and the vacuum space is at least 15 A, which 
is large enough to avoid the interaction between periodi- 
cal images. A 7x7x 1 Monkhorst-Pack grid is used for the 
sampling of the Brillouin zone during geometry optimiza- 
tion. All the atoms in the unit cell were allowed to relax, 
and the convergence of force is set to 0.01 eV/A. Addi- 
tionally, spin polarization is turned on during the relax- 
ation processes. All other calculations of accurate elec- 
tronic properties were performed using the full-potential 
linearized augmented plane-wave [21] method as imple- 
mented in the WIEN2k code [22]. It is well known that 
different XC potentials can lead, depending on the studied 
materials and properties, to results which are in very bad 
agreement with the experiment, e.g., for the bandgap of 
semiconductors and insulators which is severely underes- 
timated or even absent [23]. For this reason, the modified 
Becke-Johnson (MBJ) [24,25] potential in the framework 
of local-density approximation (LDA) [26,27] is taken to 
calculate the bandgap of a-g-B^NsC while the magnetism 
in ^-£-B 3 N 3 C is described by the PBE (the standard GGA 
for materials) potential. 

Results and discussion 

Electronic bandgap in a-g-B3N3C 

Figure 1 plots the total energy per atom against the lattice 
constant a (the lattice constant c is fixed) for g-B3N3C. 
It can be seen that the total energy of g-BsNsC as a 
function of lattice constant a has a single minimum, 
meaning that the geometrical structure would be sta- 
ble. Particularly, the charge population analysis reveals 
that the electron density around the C-N bond in a-g- 
B3N3C is much higher than that around the C-B bond 
in fi-g-BsNsC, showing that the C-N bond is relatively 
strong, which is also consistent with the experimental 
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results [1]. This strong interaction between C and N 
atoms in a-g-BsNsC directly results in a short C-N bond 
length (see the following paragraph) and can balance the 
strain of the monolayer graphite-like structure. Thereby, 
a-g-B^NsC could be a more thermodynamically stable 
topological phase against P-g-B^NsC, which has a lower 
total energy of around 0.43 eV/atom compared with 
0-S-B3N3C. 

To explore further the mechanical stability of ot-g- 
B3N3C, the optimized lattice constant a — 4.67 A is first 
obtained, as depicted in the left panel of Figure 1. Impor- 
tantly, the C-N bond length converged to 1.31 A, which 
is considerably reduced from the typical bond lengths of 
1.37 to 1.48 A in the related materials [28]. This strongly 
suggests the nature of the higher binding energy of the 
C-N bond [1], which also influences the B-N bonding 
and extends its length. The obtained B-N bond length is 
1.48 A, which is slightly bigger than the value of 1.45 A in 
h-BN. The calculated partial density of states (PDOS) is 
shown in Figure 2. The valence band is dominated by B p 
and C p states, while the conduction band is only domi- 
nated by Bp states. There, one can find the majority of Cp 
states to be semicore, lying 6 to 9 eV below the Fermi level. 
These states interact with those comprising the valence 
band with the same symmetry. As a result, there is a small 
admixture of C p and N p states close to the Fermi level. 
However, in a-g-B^NsC, one finds significant admixture 
ofCp and N p states in the semicore energy window with 6 
to 9 eV below the Fermi level. This suggests that the C p-N 
p interaction in the semicore region contribution to the C- 
N bonding is significantly more important in a-g-BzNsC 
than the C p-N p interaction close to the Fermi level. 

Now, let us look at the band structures of a-g-BzNsC, as 
given in Figure 3. It explicitly demonstrates that all three 
XC potentials give similar band structures. The highest 
occupied crystalline orbitals are located at the K point of 
the reciprocal space, while the lowest unoccupied crys- 
talline orbitals appear at the M point. This leads to an 
indirect bandgap semiconductor. To obtain the bandgap 
more accurately (see Figure 3c), the MBJ XC functional 
is used [24,25]. The bandgap of a-£-B 3 N 3 C within MBJ 
is obtained to be 1.22 eV, which nicely locates the mid- 
dle region between 0.59 and 1.80 eV for the BNC samples 
with 12.5% and 50% C contents, respectively [4]. Note that 
the C content in a-g-BsNsC is 25%. In the experiments 
[1], the absorption edges are redshifted as the C concen- 
trations increase, which shows a tunable mechanism of an 
optical bandgap in actual applications. By comparing with 
[1] where BNC with around 65% C concentration shows 
an optical bandgap of 1.62 eV, we infer that such a higher 
energy absorption edge (take into account the bandgap of 
1.22 eV in our case with 25% C concentration) arises from 
the formation of individual BN and graphene domains. In 
this way, the even distribution of C in BNC systems might 



serve as a good guide to find alternative solutions to exist- 
ing bandgap-engineered applications. Future research can 
test this prediction directly. 

In addition, the bandgaps based on the LDA and PBE 
methods are equal to 0.83 eV (see Figure 3a,b). The rela- 
tive bandgap correction from MBJ with respect to LDA, 
K = (A M bj - Alda)/A M bj with A M bj/lda being the 
bandgap, is about 32%. We can see that our calculated 1Z 
for a-g-B<$N<$C ( shows an excellent agreement with the 
16 sp semiconductors) lies within the range of 16.0% to 
100.0% (see Table 1 for details). The correction value is 
particularly very close to that in BN (25%), GaN (42%), 
A1P (37%), and A1N (25%). This finding is not surprising 
because the listed four solids have at least one element 
close to that of a-g-B^NsC in the periodic table of ele- 
ments. It means that the similar chemical circumstances 
in these materials can be well described by the same 
XC functionals. This underlyingly confirms our predic- 
tion validity and that the a-g-BsNsC might be carried out 
experimentally. 

From Figure 2, which shows the DOS of a-g-BzNzC 
we can see that the effect of the MBJ potential is to 
shift up (with respect to LDA/PBE) the unoccupied B 
2s and 2p states. Here, three major differences between 
the LDA/PBE and MBJ methods can be extracted: 

(a) Another obvious effect of MBJ potentials is to shift 
down the middle of the valence band at around —4.0 eV. 

(b) The hybridization of s and p states of dominant 
B 2s and other atoms' 2p states at the bottom of the 
valence band is very strong in MBJ calculations (denoted 
as down arrows in Figure 2). (c) The correction char- 
acter of a-g-B^NsC is much more pronounced with the 
MBJ than with the LDA/PBE, which narrows the valence 
band just below the Fermi level. We would like to stress 
that the MBJ potentials open a bandgap of 0.39 eV in 
the a-g-B^NsC model compared to the result of LDA, 
which is consistent with the orbital-dependent potential 
principle [24]. 

The optical absorption spectrum for a-g-BsNsC within 
MBJ potentials (see Figure 4c) shows a big absorption 
packet with two adjacent peaks in the range of 2.5 to 
5.0 eV, which originates from the band structure as shown 
in Figure 3c. As an indirect bandgap material, the tran- 
sition from the K point to T point is very weak as the 
momentum conservation rule is not satisfied here, and 
thus, the corresponding characteristic photoluminescence 
in the optical absorption spectrum can be negligible. This 
is also true in our case for a-g-B^NsC. The first peak cor- 
responds to the direct bandgap transition at the M point. 
The second peak comes from the larger direct gap from 
higher energy states located at the V point. More impor- 
tantly, these two peaks are not separated distinctly. This 
feature can be attributed to the even distribution of C 
atoms in a-g-B^NsC as mentioned above, which shows a 
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Figure 3 Calculated band structures for a-g-B 3 N 3 C model with the three XC potentials, (a) LDA XC potential, (b) PBE XC potential, (c) MBJ XC 
potential. 



different formation mechanism compared with the hybrid 
BNC in the experiment [1]. 

Correction effects were taken into account by adding 
a LDA correction potential in MBJ [25]. This important 



Table 1 The bandgaps and the relative bandgap correction 
71 of 1 6 sp semiconductors and the predicted a-g-BsN^C 



Solid 


LDA 


MBJ 


1Z 


Expt. 


C 


4.11 


4.93 


16.6 


5.48 


Si 


0.47 


1.17 


59.8 


1.17 


Ge 


0.00 


0.85 


100.0 


0.74 


LiF 


8.94 


12.94 


30.9 


14.20 


LiCI 


6.06 


8.64 


29.9 


9.40 


MgO 


4.70 


7.17 


34.4 


7.83 


ScN 


-0.14 


0.90 


115.6 


0.90 


SiC 


1.35 


2.28 


40.8 


2.40 


BN 


4.39 


5.85 


25.0 


6.25 


GaN 


1.63 


2.81 


42.0 


3.20 


GaAs 


0.30 


1.64 


81.7 


1.52 


AIP 


1.46 


2.32 


37.1 


2.45 


ZnS 


1.84 


3.66 


49.7 


3.91 


CdS 


0.86 


2.66 


67.7 


2.42 


AIN 


4.17 


5.55 


24.9 


6.28 


ZnO 


0.75 


2.68 


72.0 


3.44 




0.83 


1.22 


32.0 





The theoretical and experimental bandgaps (in eV) of the 1 6 sp semiconductors 
are directly taken from [25]. The relative bandgap correction 1Z, (in %) is 
calculated from the equation = (A M bj - A L da)/ Ambj, where A M bj and A L da 
are the calculated bandgaps using XC potentials MBJ and LDA. 



physical effect opens an additional bandgap by mimick- 
ing very well the behavior of orbital-dependent potentials 
and causes a rigid blueshift of the absorption spectrum 
compared with the LDA/PBE curves as shown in Figure 4. 
This explains the excellent qualitative agreement of the 
hybrid exchange-potential optical absorption spectrum 
seen in Figure 4, due to a compensation of significant 
errors within the standard DFT methods. 

Magnetism in /?-g-B3N3C 

The predicted structure of /3-JJ-B3N3C is shown in the 
right panel of Figure 1. The lattice constant of /3-/J-B3N3C 
is obtained to be 5.26 A as shown in Figure 1. The results 
show the equilibrium value cIbc = 1-52 A, which is close 
to the value of graphite-like BC3, which is 1.55 A [29]. It 
is to be noticed that all the equilibrium values, <^bn> m 
g-B^NsC are equal to 1.42 A, which is slightly less than 
the value of 1.45 A in the pristine BN sheet. This implies 
the stronger B-N bonds formed in P-g-B^NsC. Our calcu- 
lations show that the P-g-B^NsC leads to a ground state 
with a magnetic moment of 0.68 /xb. The nonmagnetic 
state is 0.07 eV higher than this ground state. 

From the band structures, we see that although both 
the pristine BN sheet and graphene are nonmagnetic, the 
P-g-B^NsC model can be spin-polarized. It is necessary 
to discuss magnetism in more detail from its electronic 
structures. Figure 5a,b,c presents the band structure and 
spin-resolved total density of state (TDOS). Remarkably, 
two bands cross the Fermi level (black and red circles 
in Figure 5a) and make the Fermi energy level occupied 
completely along the entire high-symmetry lines. In con- 
trast, the spin-down one does not possess such a strongly 
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Figure 4 The optical absorption expressed by the imaginary part of the dielectric tensor e 2 for a-g-BsNsC. The imaginary part of the 
dielectric tensor is averaged for four different directions. Three XC potentials including (a) LDA, (b) PBE, and (c) MBJ are shown in the calculations. 



localized feature, just one band (black circles) crosses the 
Fermi level monolithically as shown in Figure 5c. The 
calculated magnetic moment in fi-g-B^NsC should orig- 
inate from this asymmetric spin-dependent localization. 
The corresponding strong spin splitting can be further 
confirmed from TDOS as shown in Figure 5b. In addi- 
tion, close examination of the top valence bands (see 
Figure 6a) indicates that the strong localization mainly 
comes from the 2p atomic orbitals of C and N atoms. Here, 
we also show the Fermi surfaces of p-g-B^NzC in the first 
Brillouin zone. The unique feature of the Fermi surface 
which is almost parallel to the high-symmetry lines (see 
Figure 5d,e,f ) is the direct manifestation of the bands near 
the Fermi level in Figure 5a,c. 

It should be noticed that such magnetism is induced 
without transition metals and without external perturba- 
tions so that p-g-B^NsC behaves as the first, theoreti- 
cally predicted, metal-free magnetic material in the hybrid 
BNC system. Apparently, our finding points out a new 
direction for further related experimental investigations in 
spintronics. 

We now address the possibility of the origin of mag- 
netism in fi-g-BsNsC. Based on the orbital-resolved den- 
sity of state as shown in Figure 6a, the magnetic moment is 
mainly ascribed to the 2p orbitals of C and three N atoms. 
The spin polarization of C atom offers a magnetic moment 
of 0.18 /xb. A total magnetic moment of 0.18 /xb is shared 
equally by the 2p orbitals of three N atoms. The remain- 
ing magnetic moments distribute evenly in the interstitial 



region among N atoms. This is conceivable because fi-g- 
B3N3C, if compared with other BNC systems, has large 
interspaces between N atoms (see the 2x2 supercell in 
Figure 1 or Figure 6b). The remarkable feature is that the 
equilibrium surface density of $-g-B<$N<$Q is around 1.27 
times larger than that of a-g-B-^^Q. Owing to the large 
interspaces between atoms, the hydrogen storage in /3- 
g-B^NsC may be expected [30]. The detailed description 
of hydrogen storage related to P-g-B^NsC is beyond the 
scope of this work. 

Concerning the special carbon atom in p-g-BsNsC the 
counting of four valence electrons is as follows: Three 
electrons participate in the sp 2 hybrid orbital, which forms 
a planar structure. The remaining one electron is then 
redistributed in the whole unit cell due to the enhanced B- 
N covalent bond (with shorter bond length compared with 
the value in pristine h-BN), which makes the magnetic 
properties more complicated. From the fourth electron, 
only 18% of the electron still fills the n -orbital of the C 
atom and contributes a magnetic moment of 0.18 /xb. This 
is in fairly good agreement with theoretical description 
of the 7i -orbital state to the local magnetic moment of 
approximately 0.3 /xb in the graphene-like systems [31,32]. 
Around 32% of the fourth electron mainly resides at the 
interstitial region, which plays a crucial role in fi-g-B-sN-sC 
as follows: (a) enhances the B-N covalent bond, (b) pro- 
vides the main interstitial magnetic moment of 0.32 /xb, 
and (c) promotes the 2p z of the N atom spin-polarized 
slightly with a magnetic moment of 0.06 /xb per N atom. 
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The remaining percentage of the fourth electron acts as 
the conduction electrons and makes the system metallic, 
which dominates the mechanism of ferromagnetic order- 
ing in P-g-B^NsC. In the case of P-g-B^NsC we can see 
that the electron spin at the localized n -orbital state of 
C and N atoms as well as the interstitial region compels 
two energy bands localized strictly along the entire high- 
symmetry lines, i.e., the T-M-K-r line (see Figure 5a). 
Thus, the RKKY interaction [10,33] among the magnetic 



sites through the residual conduction electrons forms a 
spin ordering in these orbitals, which is the physical origin 
of ferromagnetism in P-g-B^-^C, Figure 6b,c respectively 
plot the top and side views of the 3D iso-surfaces for 
net magnetic charge density in the xy plane. This find- 
ing is insightful, and three major points deserve comment: 
(a) The C site is more spin-polarized as compared with 
each N site, (b) The dumbbell-like magnetic moment dis- 
tribution along the z direction implies that the 2p z orbital 




E-E F (eV) 

Figure 6 PDOS and spin density for ^-g-B 3 N 3 C. (a) PDOS on the interstitial region and on the 2s and 2p orbitals of C and N atoms in /3-g-B3N 3 C 
(b) The 3D iso-surface plot of spin density for the 2x2 supercell at the value of 0.04 e/A 3 . (c) A side view of spin density corresponding to (b). 
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becomes partially filled with one spin-up electron, (c) The 
induced moments are ferromagnetic coupled between the 
N and C sites based on the RKKY exchange interaction 
model as mentioned above. 

Conclusions 

In summary, we have predicted a novel crystalline mate- 
rial g-B^NsC, which displays two distinct electronic prop- 
erties where the selective bonding type of the C atom is a 
key parameter for future industrial processes. a-g-BzN^C 
is a semiconductor, while fi-g-B^NsC behaves like a metal 
and holds a magnetic moment of 0.68 /xb. Importantly, 
compared with the hybrid BNC, g-I>3N3C is proposed to 
have a simple structure, which can be applied in various 
fields due to its unique properties. 
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